Quark mass dependence of p and a Mesons from Dispersion Relations 

and Chiral Perturbation Theory 
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We use the one-loop Chiral Perturbation Theory Trvr-scattering amplitude and dispersion theory 
in the form of the inverse amplitude method to study the quark mass dependence of the two lightest 
resonances of the strong interactions, the /o(600) (a) and the p meson. As main results we find 
that the pmr coupling constant is almost quark mass independent and that the p mass shows a 
smooth quark-mass dependence while that of the a shows a strong nonanalyticity. These findings 
are important for studies of the meson spectrum on the lattice. 
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Although studied theoretically as well as experimen- 
tally for many years, the spectrum of the lightest reso- 
nances in QCD is still not understood from first princi- 
ples. The only known way to extract nonperturbative 
quantities from QCD is the use of lattice QCD. However, 
current calculations are typically still done for relatively 
high quark masses (see, e.g., Refs. [11,0). Thus, in order 
to make contact with experiment, appropriate extrapola- 
tion formulas need to be derived. This is typically done 
by using chiral perturbation theory (ChPT), the low en- 
ergy effective theory of QCD 0, ChPT predictions 
are model independent and, in particular, provide, as an 
expansion, the dependence of observables on the quark 
masses (or equivalently the pion mass). The aim of this 
Letter is to predict the quark-mass dependence of the a 
and the p mesons from basic principles, namely, using 
ChPT to next-to-leading order (NLO), unitarity, and 
analyticity in the form of dispersion theory using the in- 
verse amplitude method (JAM) [5|]. It is obtained from a 
subtracted dispersion relation of the inverse amplitude, 
whose imaginary part in the elastic region is known ex- 
actly from unitarity. All dependences on QCD parame- 
ters appear through the ChPT expansion, which is used 
to calculate the low energy subtraction points and the 
left cut. Hence, up to a given order in ChPT, the ap- 
proach has no model dependences. For this work, we will 
use the NLO SU{2) elastic JAM whose parameters are 
fitted to TTvr scattering data, and therefore, our results 
are model independent only up to NLO ChPT. 

The use of dispersive methods also allows for a 
straightforward extension to the second Ricmann sheet of 
the complex plane where poles associated to resonances 
occur. In this way both a and p poles appear naturally 
without any further assumptions and have the correct 
dependence on QCD parameters up to the order of the 
ChPT expansion used in the lAM. For instance, from a 
study of the leading 1 /Nc behavior of the amplitudes it 
was possible to conclude that the p is mostly of qq nature 
whereas the cr is predominantly non-gq f^. 

In this Letter, we study how this different structure 
encoded in the ChPT parameters gets reflected in the 
quark mass dependence. As we will see, also with respect 



to the quark-mass dependence, the p behaves differently 
than the a. We will show that for sufficiently large quark 
(pion) masses, both states become stable poles on the 
physical sheet; however, this limit is approached very 
differently, only in part due to their different quantum 
numbers. 

Another motivation for this study is the Anthropic 
Principle i-e., the need for a subtle fine tuning of 
various parameters of the Standard Model. In order to 
allow for an efficient triple a process, necessary for the 
production of carbon, a NN interaction within 2 % of the 
known strength is necessary. Since the a meson plays a 
central role in the NN interaction, this leads to bounds 
for the sigma mass and, correspondingly, for the quark 
masses 0. This issue is particularly exciting, if the fun- 
damental constants were time dependent, as claimed in 
Ref. 0. In this context the quark mass dependence of 
the cr was studied in Ref. p^ . 

We are interested in tttt elastic amplitudes projected on 
partial waves tjj of definite isospin / and total angular 
momentum J. For simplicity, we will drop the IJ labels. 
Note that the a (p) resonance appears as a pole in the 
second Riemann sheet of the partial wave (/, J) — (0, 0) 
((1, 1)). Elastic unitarity implies: 



luit{s) = a{s)\t{s)\^' 



Im [t(s) 



-a(.), (1) 



where s is the usual Mandelstamm variable for the total 
energy, a(s) = p/{2y/s) and p is the center of mass mo- 
mentum. Consequently, the imaginary part of the elastic 
inverse amplitude is known exactly. 

In this Letter, we focus on the two lightest resonances 
of QCD, the p(770) and the /o(600). It is therefore 
enough to work with the two lightest quark flavors u, d 
in the isospin limit of an equal mass m = (m„ -I- md)/_2. 
The pion mass is given by an expansion m J ^ m + ... . 
Therefore, studying the quark mass dependence is equiv- 
alent to studying the ttItt dependence. 

ChPT amplitudes are obtained as a series expansion 
t = t2 + t4 + ... with tk = 0{p^), where p stands gener- 
ically either for pion momenta or masses. The leading 
order, ^2, which is obtained at tree level from the O(p^) 
Lagrangian, is just a polynomial fixed by chiral symme- 
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try in terms of the pion mass and its decay constant /tt. 
The NLO term ^4 has both one-loop contributions from 
the 0{p^) Lagrangian and tree level contributions from 
the 0{p^) Lagrangian. The latter depend on a set of 
low energy constants (LEG), denoted by that absorb 
the one-loop divergences through renormalization. Since 
these constants are the coefficients of the energy and mass 
expansion the li have no quark mass dependence. In con- 
trast, /tt gets renormalized at NLO and thus depends 
explicitly on the pion mass. 

The ChPT series, being an expansion, satisfies unitar- 
ity, Eq.([T]), just perturbatively: 



{s — SA)/iz — SA) — {s~ S2)/{z~ S2)- Altogether, we find 



Imt4(s) = CT{s)\t2{sf 



Im 



U{s) 



= a{s), (2) 



and cannot generate poles. Therefore the resonance re- 
gion lies beyond the reach of standard ChPT. However, it 
can be reached by combinin g C hPT with dispersion the- 
ory either for the amplitu de llll| or the inverse amplitude 
through the lAM |l3,[ll,[l3~^ 

3. The lAM uses the ChPT series to generate reso- 
nances in meson-meson scattering. In the elastic case, 
the lAM follows from dispersion theory [l^, due to 
the fact that t and 1/t have an almost identical analytic 
structure: both have a "physical cut" from threshold to 
00 and a "left cut" from —00 to s = 0. In addition, the in- 
verse amplitude may have a pole whenever the amplitude 
vanishes. For the p-channel, this happens at threshold 
to all orders, but for scalar waves it occurs at the so- 
called Adler zero sa, that lies on the real axis below 
threshold, thus within the ChPT region of applicability. 
Its position can be obtained from the ChPT series, i.e., 
Sa = S2 + S4 + where t2 vanishes at S2, ^2 + ^4 at 
S2 + S4, etc... Thus we can write a dispersion relation for 
the inverse amplitude. 



1 

t{sj 



SA 



dz ■ 



Iml/<(z) 



RC 



[z - sa){z - s) 



LCi+PCi, 



(3) 

where "LC" stands for a similar integral over the left 
cut. To ensure convergence we have made a subtraction 
precisely at sa- Since ^2 is real on the real axis, we can 
similarly write 



his)_ 
Hsf 



S2 



dzp^^^}M^+LCu+PC... 
RC {z-S2){z-s) 17 ^ 

(4) 

We can now use unitarity, Eqs.([T|) and to find that 
the integrand numerators are exactly opposite. Since 
the LC integral is weighted at low energies, using NLO 
ChPT, we also find LC[l/t) ~ -LC{ti/tl). Conse- 
quently, crossing symmetry is satisfied just up to NLO. 
In the above relations, PCi/t and PC^^/f^ stand for the 
contributions of a double and triple pole respectively, 
which can be easily calculated within ChPT. Their con- 
tribution is only relevant around the Adler zero, where 
they diverge. Finally, we can use ChPT to approximate 



t2{s)-t4{s)+A^^^^^f{s) 

, , . {s2-Sa){s ~S2) [t'2{s2)-t'i{s2)] 
ti{S2) 



(5) 



~SA 



The usual lAM is recovered for ^ which holds 

exactly for all partial waves but the scalar ones. In the 
original lAM dispersive derivation [H, E^l A"^^^^'^ was 
neglected, since it formally yields a NNLO contribution. 
However, due to neglecting A^^^^ , the lAM has a spu- 
rious pole close to its Adler zero, which is only correct 
to LO ChPT. As we will see, for large m^r the a pole 
splits into two virtual poles below threshold, one of them 
moving towards zero and eventually approaching the spu- 
rious pole of the lAM. Thus, although the Adler zero 
(~ rn^/-\/2) is very deep in the subthreshold region, we 
used the modified lAM (mlAM) [1, which has no 
spurious pole, and reproduces the Adler zero up to NLO. 
Switching from the lAM to the mlAM influences only the 
mentioned second a pole, and only when it is very close 
to the spurious pole (this occurs when A/^ < 1.5to^). 
Besides this, the I AM and ml AM results are essentially 
the same. For subtractions made at different low en- 
ergy points, ^"i-f-^A/ acquires additional NNLO terms, 
but their effect is negligible 

Therefore, up to a given order in ChPT, the elastic 
(m)IAM is built in a model independent way from the 
first principles of unitarity and analyticity in the form of 
a dispersion relation. The ChPT series is used only on 
the Adler zeros and the left hand cut which is heavily 
weighted at low energies, thus well within its region of 
applicability. It is a dispersion integral for the inverse 
amplitude that allows us to study the resonance region. 

Although remarkably simple, the (m)IAM amplitudes 
satisfy elastic unitarity, Eq.([T]), exactly, and provide a 
very good description of meson-meson scattering data 
simultaneousl y i n the resonance and low energy re- 
gions [H, [HjHI]- Furthermore, the lAM generates the 
poles in the second Riemann sheet associated to the res- 
onances, namely a and p. This description is obtained 
with values of the LEC compatible with those of standard 
ChPT Actually, the ChPT series up to NLO is 

recovered when Eq. ([5]) is reexpanded at low energies. 
Thus, after a fit to data we can modify to^ and follow 
the poles associated to p and a on the second sheet. 

As long as they fall within their uncertainties, the pre- 
cise values of the LECs Z3 and l\ are not very relevant for 
this study. We take from |4|] \Q% = 0.8 ± 3.8, IQ% = 
6.2 ± 5.7. Then we fit the mlAM to data up to the 
resonance region and find 10^Z[ = —3.7 ± 0.2, 10'^Z2 = 
5.0 ± 0.4. All these LEC are evaluated at = 0.77 GeV. 

The values of considered should fall within the 
ChPT range of applicability and allow for some elastic tttt 
regime to exist below KK threshold. Both criteria are 
satisfied, if < 0.5 GeV, since we know SU(3) ChPT 
still works fairly well with such a kaon mass, and because 
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FIG. 1: Movement of the cr (dashed lines) and p (dotted lines) 
poles for increasing pion masses (direction indicated by the 
arrows) on the second sheet. The filled (open) boxes denote 
the pole positions for the a (p) at pion masses m-n — 1, 2, and 
3 X m^^^^, respectively. Note, for nirr — 3m^^^^ three poles 
accumulate in the plot very near the nn threshold. 




for m^r — 0.5 GeV, the kaon mass becomes ~ 0.6 GeV, 
leaving a 0.2 GeV gap to the two-kaon threshold. For 
larger values of a coupled-channel lAM is needed, 
which is feasible, but lies beyond our present scope, and 
lacks a dispersive derivation. 

Fig. [T] shows, in the second Riemman sheet, the p and 
a poles for the physical m^, and how they move as 
increases. Note that, associated to each resonance, there 
are two conjugate poles that move symmetrically on each 
side of the real axis. In order to see more clearly that 
all poles move closer to the two-pion threshold, which is 
also increasing, all quantities are given in units of 
so that the two-pion threshold is fixed at ^/s = 2. Let 
us recall that, for narrow resonances, their mass Af and 
width F are related to the pole position in the lower half 
plane as y/Spoie — M — iT/2 and customarily this no- 
tation is also kept for broader resonances. Hence, both 
Fcr and Fp decrease for increasing m^. In particular, Fp 
vanishes exactly at threshold where one pole jumps into 
the first sheet, thus becoming a traditional stable state, 
while its partner remains on the second sheet practically 
at the very same position as the one in the first. In 
contrast, when Mo- reaches the two-pion threshold, its 
poles remain on the second sheet with a non-zero imagi- 
nary part before they meet on the real axis and become 
virtual states. As TOtt increases further, one of those vir- 
tual states moves towards threshold and jumps onto the 
first sheet, whereas the other one remains in the second 
sheet. Such an analytic structure, with two very asym- 
metric poles in different sheets of an angular momentum 
zero partial wave, is a strong indication for a prominent 
molecular component [l^, [l3|. Differences between P- 
wave and S-wave pole movements were also found within 
quark models [1^, the latter also showing two second 
sheet poles on the real axis below threshold. 

In the upper panel of Fig. [2] we show the depen- 



FIG. 2: m-^ dependence of resonance masses (upper panel) 
and widths (lower panel) in units of the physical values. In 
both panels the dark (light) band shows the results for the a 
(p). The width of the bands reflects the uncertainties induced 
from the uncertainties in the LEG. The dotted line shows 
the a mass dependence estimated in Ref. Q|. The dashed 
(continuous) line shows the dependence of the a (p) width 
from the change of phase space only, assuming a constant 
coupling of the resonance to tttt. 



dence of and Mp normalized to their physical values. 
The bands cover the LEG uncertainties. Note, that sig- 
nificant, additional uncertainties may emerge at the two 
loop level for pion masses larger than 0.3 GeV — see, 
e.g., Ref. [l^. We see that both masses grow with in- 
creasing m,r, but the rise of is stronger than that of 
Mp, and again we see that around m^r — 0.33 GeV the a 
state splits into two virtual states with different behav- 
ior. The upper branch moves closer to threshold and thus 
has the biggest influence in the physical region, eventu- 
ally jumping to the first Riemann sheet. Note that the 
m^r dependence of is much softer than that suggested 
in the model of Q, shown as the dotted line, which in 
addition does not show the virtual pole splitting. 

In the lower panel of Fig. [2] we show the rriTr depen- 
dence of Fg- and Fp normalized to their physical values. 
The decrease in Fp is largely kinematical, following re- 
markably well the expected reduction from phase space 
as m^r and Mp increase. In other words, the effective cou- 
pling of the p to TTTT is almost m,r independent. This was 
assumed in the analysis of Ref. [20|; however, so far this 
assumption has not been supported by theory. In sharp 
contrast to this behavior is the one of Fg-. This suggests 
a strong pion mass dependence of the a coupling to two 
pions, necessarily present for molecular states [l7l.l2ll|. 

Fig. [3] is a comparison of our results for the m-^ de- 
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FIG. 3: The grey band shows the m,r dependence of 
mass from the lAM versus recent lattice resuhs from UJ 
dashed hne is the lAM result for Nc = 10. 
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pendence of Mp with some recent unquenched lattice re- 
sults [l| , which deserves several words of caution. In par- 
ticular, our approach only ensures the to^ dependence 
contained in the NLO ChPT series — it is, e.g., lack- 
ing terms of order m^p^ and higher. In addition, Mp is 
the 'pole mass' which, particularly for physical values, is 
deep in the complex plane, while, due to the finite lat- 
tice volume, the minimum energy with which pions are 
produced on the lattice is larger than the resulting Mp. 
In our formalism, we can mimic a narrow p by increas- 
ing the number of colors We also show the result 
of rescaling the lAM ChPT amplitudes from Nc = i 
to Nc = 10, which effectively reduces Tp by a factor of 
3/10. Although this narrowing effect is not exactly the 
same as that on the lattice - it is more like quenching the 
lattice results, since the large expansion actually sup- 



presses quark loops - it is encouraging that making the 
p artificially narrower yields a better agreement with the 
quark-mass dependence of the lattice data. With these 
caveats in mind our results are in qualitative agreement 
with the lattice results. Following Ref. one may write 
Mp — Mp + ami + 0(mj) , where the a parameters 
are expected to be of order one and 0.65 GeV < M^ < 
0.80 GeV. This is confirmed by a fit to lattice data [l|. 
From our approach we predict M° = 0.735 ±0.0017 GeV. 
Furthermore, the lAM reproduces Mp at the physical 
value of mTT, where higher orders in the p mass formula 
are ~ 15% [22j. Within that uncertainty, we thus get 
the prediction ci = 0.90 ± 0.11 ± 0.13 GeV"\ Although 
the quark-mass dependence of our calculation is steeper 
than that of Ref. [l| , the extracted values are still con- 
sistent with the expectations mentioned above — again 
we remind the reader that the rriTr dependence included 
is correct only to NLO in ChPT. 

To summarize, we presented a prediction for the quark 
mass dependence of the lightest resonances in QCD, 
namely the p and the a meson based on chiral pertur- 
bation theory at next -to-leading order together with the 
inverse amplitude method. We showed that the mass of 
the p has a very smooth m^r dependence and its coupling 
to TTTT is almost quark-mass independent — in Ref. [20| . 
this was only assumed without further evidence. The 
mass of the ct, on the other hand, shows a pronounced 
non-analyticity when is varied. In addition, its ef- 
fective coupling to TTTT is strongly TOjt dependent. This is 
interpreted as additional evidence for a significant molec- 
ular admixture in the sigma, consistent with previous 
analyses, and will be important for future chiral extrap- 
olations of lattice data for s-wave resonances. 

We thank G. Colangelo, S. Diirr, V.V. Flambaum, U.- 
G. MeiBner and A. Rusetsky for their useful comments. 
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